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Abstract: An ecologically industrial system is such an industrial system in harmony with 
its environment, especially the natural environment. The main purpose of this paper is to 
show how to establish a mathematical model for such systems by combinatorics, and find 
its topological characteristics, which are useful in industrial ecology and the environment 
protection. 
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§1. Introduction 

Usually, the entirely life cycle of a product consists of mining, smelting, production, storage, 
transporting, use and then finally go to the waste, • • • , etc.. In this process, a lot of waste gas, 
water or solid waste are produced. Such as those shown in Fig.l for a producing cell following. 
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Fig.l 

In old times, these wastes produced in industry are directly discarded to the nature without 
disposal, which brings about an serious problem to human beings, i.e. , environment pollution 
and harmful to our survival. For minimizing the effects of these waste to our survival, the growth 
of industry should be in coordinated with the nature and the 3 R rule: reduces its amounts, reuses 
it and furthermore, into recycling, i.e., use these waste into produce again after disposal, or let 
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them be the materials of other products and then reduce the total amounts of waste to our life 
environment. An ecologically industrial system is such a system consisting of industrial cells in 
accordance with the 3 R rule by setting up one or more waste disposal centers. Such a system 
is opened. Certainly, it can be transferred to a closed one by letting the environment as an 
additional cell. For example, series produces such as those shown in Fig. 2 following. 




product3 



Fig-2 

Generally, we can assume that there are Pi, P 2 , ■ ■ ■ , P m products (including by-products) 
and W[ , W' 2 , • • • , W s wastes after a produce process. Some of them will be used, and some will 
be the materials of another produce process. In view of cyclic economy, such an ecologically 
industrial system is nothing else but a Smarandachely multi-system. Furthermore, it is a 
combinatorial system defined following. 

Definition 1 . 1 ( [1] , [2] and [9]) A rule in a mathematical system (E;P) is said to be Smaran- 
dachely denied if it behaves in at least two different ways within the same set E, i.e., validated 
and invalided, or only invalided but in multiple distinct ways. 

A Smarandachely system (E;P) is a mathematical system which has at least one Smaran- 
dachely denied rule in P. 

Definition I.2([I],[2] and [9]) For an integer m > 2, let (Si; Pi), (E 2 ;P 2 ), •••, (S m ;P m ) 
be m mathematical systems different two by two. A Smarandache multi-space is a pair (S ;1Z) 
with 

7n m 

E = |^J E, , and 1Z = |^J P,; . 

2=1 i—1 

Definition I.3([I],[2] and [9]) A combinatorial system is a union of mathematical systems 
(Ei;Pi),(E 2 ;P 2 ), •••, (E m ;P m )/or an integer m, i.e., 

m m 

*b = ( U s*; U ftO 

2 = 1 2=1 

with an underlying connected graph structure G, where 



V(G) = {Ei, E 2 , • • • ,E m }, 

E(G) = { (Ej, Ej) | Ej p| Ej ^ 0, 1 < i,j < to}. 
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The main purpose of this paper is to show how to establish a mathematical model for 
such systems by combinatorics, and find its topological characteristics with label equations. In 
fact, such a system of equations is non-solvable. As we discussed in references [3]- [8], such a 
non-solvable system of equations has significance also for things in our world and its global 
behavior can be handed by its G-solutions, where G is a topological graph inherited by this 
non-solvable system. 

§2. A Generalization of Input-Output Analysis 

The 3R rule on an ecologically industrial system implies that such a system is optimal both in 
its economical and environmental results. 

2.1 An Input-Output Model 

The input-output model is a linear model in macro-economic analysis, established by a economist 
Leontief as follows, who won the Nobel economic prize in 1973. 

Assume these are n departments Di,D 2 ,--- i D n in a macro-economic system A? satisfy 
conditions following: 

(1) The total output value of department Dj is a Among them, there are Xij output 
values for the department Dj and di for the social demand, such as those shown in Fig.l. 

(2) A unit output value of department Dj consumes Uj input values coming from depart- 
ment Di. Such numbers iy, 1 < i,j < n are called consuming coefficients. 




Fig. 2 

Therefore, such an overall balance macro-economic system ££ satisfies n linear equations 

n 

Xi — ^ ^ Xij T di 

3 = 1 



(1) 
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for integers 1 < i < n. Furthermore, substitute tij = Xij/xj into equation (10-1), we get that 

n 

Xj — ^ ^ tij Xj T di (2) 

i= i 

for any integer i. Let T = [tij\ nxn , A = I nxn - T. Then 

Ax = d, (3) 

from (2), where x = (xi, X 2 , • • • , x n ) T , d = ( d\ , c ?2 , ■ • • , d n ) T are the output vector or demand 
vectors, respectively. 

For example, let ££ consists of 3 departments D\, D 2 , D 3 , where I?i=agriculture, D 2 = 
manufacture industry, Z? 3 =service with an input-output data in Table 1. 



Department 


Dr 


d 2 


d 3 


Social demand 


Total value 


Di 


15 


20 


30 


35 


100 


d 2 


30 


10 


45 


115 


200 


d 3 


20 


60 


/ 


70 


150 



Table 1 

This table can be turned to a consuming coefficient table by tij = Xij jxj following. 



Department 


Di 


d 2 


d 3 


D x 


0.15 


0.10 


0.20 


d 2 


0.30 


0.05 


0.30 


d 3 


0.20 


0.30 


0.00 



Thus 



Table 2 





’ 0.15 


0.10 


0.20 ’ 




0.85 


-0.10 


-0.20 ’ 


T = 


0.30 


0.05 


0.30 


i A — d 3x3 — T — 


-0.30 


0.95 


-0.30 




_ 0.20 


0.30 


0.00 _ 




_ -0.20 


-0.30 


1.00 



and the input-output equation system is 



! 0.85xi — 0.10x2 — 0.20x’3 = d\ 

— 0.30xi + 0.95x2 — 0.30x3 = d 2 
— 0.20xi — 0.30x2 + x — 3 = g ?3 



Solving this linear system of equations enables one to find the input and output data for 
economy management. 
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2.2 A Generalized Input-Output Model 

Notice that our WORLD is not linear in general, i.e. , the assumption tij = Xij/xj does not hold 
in general. A non-linear input-output model is shown in Fig. 3, where x = (xu,X 2 i,- • • ,x n i), 
Di,D 2 ,--- ,D n are n departments, SD=social demand. Usually, the function F(x) is called 
the producing function. 




Fig. 3 



In this case, an overall balance input-output model is characterized by equations 

n 

Fi(x) =^2xij + di (4) 

i= i 



for integers 1 < * < n, where Fjf'x) may be linear or non-linear and determined by a system of 
equations such as those of ordinary differential equations 



1 < i < n 



Fj ^ + a\Fj ) + ••• + a n ~\Fi + a n — 0 



Fi| t= o — <Po i Fi 



(i) 



it=o 



= <Pi, 






it=o 



— (fin— 1 



(< OES n ) 



or 



1 < i < n 



ori - u u „ 

o , -F* 1 *^1 5 i %n— 1 1 Pi 1 

Ot 

Fi\t=t 0 = ipo(xi,X2,- ■ ■ ,X n -l) 



,Pn- l) 



(PES 1 ) 



which can be solved by classical mathematics. However, the input-output model with its gener- 
alized only consider the consuming and producing, neglected the waste and its affection to our 
environment. So it can be not immediately applied to ecologically industrial systems. However, 
we can generalize such a system for this objective by introducing environment factors, which 
are discussed in the next section. 
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$3. A Topological Model for Ecologically Industrial Systems 



The essence of input-output model is in the output is equal to the input, i.e. , a simple case of the 
law of conservation of mass: a matter can be changed from one form into another, mixtures 

can be separated or made, and pure substances can be decomposed, but the total amount of 
mass remains constant. Applying this law, it needs the environment as an additional cell for 
ecologically industrial systems and replaces the departments Di, 1 < i < n by input materials 
Mi, 1 < i < n or products Pk, 1 < k < m, and SD by Wj, 1 < i < s = wastes, such as those 
shown in Fig. 4 following. 




Fig. 4 

In this case, the balance input-output model is characterized by equations 

n s 

F i (x) = H W * ( 5 ) 

j = 1 2 = 1 

for integers 1 < i < n. We construct a topological graphs following. 

Construction 3.1 Let ^ (t) be an ecologically industrial system consisting of cells C\(t), C^it), 
■ ■■ ,Ci(t), R the environment of ^ . Define a topological graph G[ c fl\ of ^ following: 

V(G[J ? ]) = ,Ci{t),R}; 

E{G[/]) = {( Ci(t),Cj(t )) if there is an input from Cjft) to Cj(t), 1 < i,j < Z} 

( J{ (Ci (t) , R) if there are wastes from Ci(t ) to R, 1 <i< l}. 

Clearly, G\J?\ is an inherited graph for an ecologically industrial system J? . By the 
3R rule, any producing process X n of an ecologically industrial system is on a directed cycle 
C = (A, :i , X i2 , • • • , Xi k ), where X l} G {Ci, 1 < j < l\R}, such as those shown in Fig. 5. 
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XiAt) 




O * o 



x h {t) • • • x ih {t) 



Fig. 5 

Such structure of cycles naturally determined the topological structure of an ecologically indus- 
trial system following. 

Theorem 3.2 Let ^ (t) be an ecologically industrial system consisting of producing cells 
C 2 A), • • • underlying a graph G[ c /(t)\. Then there is a cycle- decomposition 

t 

G[j?(t)] = \JC ki 

i= 1 

for the directed graph G [jP{t)\ such that each producing process Gift), 1 < i < l is on a directed 
circuit C ki for an integer 1 < i < t. Particularly, G[^(t)\ is 2-edge connectness. 

Proof By definition, each producing process Ci(t) is on a directed cycle, which enables us 
to get a cycle-decomposition 

t 

GG[j?(t)] = \JC ki . □ 

2 = 1 

Thus, any ecologically industrial system underlying a topological 2-edge connect graph 
with vertices consisting of these producing process. Whence, we can always call G-system for 
an ecologically industrial system. Clearly, the global effects of Gi-system and G 2 -system are 
different if G\ qk. G 2 by definition. Certainly, we can also characterize these G-systems with 
graphs by equations (5) following. 

Theorem 3.3 Let consisting of producing cells Gift), C 2 A), ■ ■ ■ ,Ci(t ) underlying a graph 
G[/{t)}. Then 



F v (x uv ,uG N~ [J{t)] (v))= (-l) 5(v ’ w) x vw 

weN oumG) 

with 8(v,w) = 1 if x vw =product, and —1 if x vw =waste, where are ^ ie or 

our-neighborhoods of vertex v in G[^ (<)]. 

Notice that the system of equations in Theorem 3.3 is non-solvable in R A+1 with A the 
maximum valency of vertices in G (t)\. However, we can also find its G \jA (f)] -solution in 
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R a+1 (See [4]- [6] for details), which can be also applied for holding the global behavior of such 
G-systems. Such a G\^{t )]- solution is not constant for Ve € E(G [^ (£)])• For example, let a 
G-system with G=circuit be shown in Fig. 4. 




G 6 



Fig. 5 

Then there are no wastes to environment with equations 

F v (x v J = x Vi+1 , 1 < i < 6, where i mod6, i.e. , 

F Vi F Vi+1 ■ ■ ■ F Vi+6 = 1 for any integer 1 < i < 6. 

If F Vi is given, then solutions x Vi , 1 < * < 6 dependent on an initial value, for example, 
x V \\t=o , i- e -> one nee ds the choice criterions for determining the initial values x Vi \ t _ (y Notice 
that an industrial system should harmonizes with its environment. The only criterion for its 
choice must be 

optimal in economy with minimum affection to the environment, or approximately, maxi- 
mum output with minimum input. 

According to this criterion, there are 2 types of G-systems approximating to an ecologically 
industrial system: 

(1) Optimal in economy with all inputs (wastes) W ri , W r2 ,- ■ ■ , W Ts licenced to R; 

(2) Minimal wastes to the environment, i.e., minimal used materials but supporting the 
survival of human beings. 

For a G-system, let 

c~ = ^2 c(x uv ) and c+ = ^ {-^) 5(v ' w) c{x vw ) 

be respectively the producing costs and product income at vertex v £ V(G). Then the optimal 
function is 
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A (G) = £ (4-c-) 

vev(G) 



( 



= Y E (-^) 5{v ’ w) c(x vw )~ y ■ 



veV(G) \t«6JV+ [</(t)] (t>) uiN GlJfm] (v 

Then, a G-system of Types 1 is a mathematical programming 

Y A (G) but ^ W H < W-; c 



max 



rU 
' ri > 



vev(G) vev(G) 

where is the permitted value for waste W r i to the nature for integers 1 < i < s. Similarly, 
a G-system of Types 2 is a mathematical programming 

min Y^ W r i but all prodcuts P > P L , 



ve v(G) 

where P L is the minimum needs of product P in an area or a country. Particularly, if W^[ = 0, 
i.e., an ecologically industrial system, such a system can be also characterized by a non-solvable 
system of equations 

Fv (%UV ? ^ ^ N G[Jf(t)]( V )) - Y X V yj for \/v G V(G). 

WeN GlJ?(. t)]G) 
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